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ABSTRACT 


An extension of the two-dimensional formulation developed last year is presented 
for a three dimensional body of revolution. With the introduction of a Fourier expansion of 
the vector electric and magnetic fields, a coupled two-dimenional system is generated and 
solved via the finite-element method. As before, an exact boundary condition is employed 
to terminate the mesh and the FFT is used to evaluate the boundary integrals for low O(n) 
memory demand when an iterative solution algorithm is used. Again, by virtue of the finite 
element method the algorithm is applicable to structures of arbitrary material composition. 

Several improvements to our two-dimensional algorithm are also described. These 
include (1) modifications for terminating the sh at circular boundaries without distorting 
the convolutionality of the boundary -.rr.egr-.s, (2) the development of our own non- 
proprietery mesh generation routines for two-dimensional applications, (3) the development 
of preprocessors for interfacing SDRC IDEAS@ with the main algorithm, and (4) the 
development of post-processing algorithms based on the public domain package GRAFIC 
to generate 2D and 3D gray level and color field maps. 


OBJECTIVE 

The objective of this task is to develop innovative techniques and related software 
for scattering by three dimensional composite structures. The proposed analysis is a hybrid 
finite element-boundary integral method formulated to have an O(n) memory demand. This 
low storage is achieved by employing the FFT to evaluate all boundary integrals and 
resorting to an iterative solution algorithm. Particular emphasis in this task is the 
generation of software applicable to airborne vehicles and the validation of these by 
comparison with measured and other reference data. Because the approach is new, a step 
by step development procedure has been proposed over a three-year period. During the 
first year the technique was developed and implemented for two-dimensional composite 
structures. Support software for the two-dimenional analysis such as pre- and post- 
processor routines were developed during the second year and a formulation was also 
developed and implemented for three-dimensional bodies of revolution. Finally, during the 
third year, we will develop, implement, and test the method for arbitrary three dimensional 
structures. 
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BACKGROUND 


Interest in three-dimensional (3-D) methods has increased in recent years, however, 
the associated demands in computation time and storage are often prohibitive for electrically 
large 3-D bodies. Vector and concurrent (i.e. hypercube, connection, etc.) computers are 
beginning to alleviate the first of these demands, but a minimization of the storage 
requirements is essential for treating large structures. 

The traditional Conjugate Gradient Fast Fourier Transform (CGFFT) method [ 1 ]- 
[4] is one such frequency domain solution approach which requires O(n) storage for the 
solution of n equations. This method involves the use of FFTs whose dimension equals 
that of the structure under consideration [5]-[7] and, therefore, demands excessive 
computation time when used in an iterative algorithm. Also, the standard CGFFT requires 
uniform rectangular gridding that unnecessarily includes the impenetrable portions of the 
scatterer. With these issues in mind, a new solution approach is propsed for solving 
scattering problems. The proposed method will be referred to as the Finite Element- 
Conjugate Gradient Fast Fourier Transform (FE-CGFFT) method. 

During last year's effort the FE-CGFFT method was developed for two- 
dimensional scatterers where the finite element mesh was terminated at a rectangular box. 
Inside the box boundaries, Helmholtz equation is solved via the finite element method and 
the boundary constraint is obtained by an appropriate integral equation which implicitly 
satisfies the radiation condition. Along the parallel sides of the box, this integral becomes a 
convolution and is, therefore, amenable to evaluation via the FFT. The dimension of the 
required FFT in this hybrid method is one less than the dimensionality of the structure thus, 
making it attractive for 3-D simulations. Also, because it incorporates the finite element 
method, the FE-CGFFT formulation remains valid regardless of the structure's geometry 
and material composition. 

The proposed method described in the University of Michigan Report 025921-6-T 
(see also [8]) is similar to the moment method version developed by Jin [9]. Jin's method 
was in turn based on work published in the early 70 s by McDonald and Wexler [10] who 
introduced an approach to solve unbounded field problems. The proposed method is also 
similar to other methods (a few of which will be mentioned here), neither of which 
provides a storage reduction comparable to the proposed FE-CGFFT method. The 
unimoment method [11] uses finite elements inside a fictitious circular boundary and an 
eigenfunction expansion to represent the field in the external region. The coefficients of the 
expansion are then determined by enforcing field continuity at the finite element (FE) mesh 
boundary. The coupled finite element-boundary element method [12] uses the finite 
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element method within the boundary and the boundary element method to provide the 
additional constraint at the termination of the mesh. Unlike the proposed method, the 
solution in [12] was accomplished by direct matrix inversion (as in [9]), and the outer mesh 
boundary is not rectangular to take advantage of the FFT for the evaluation of the boundary 
integrals. 


PROGRESS 

The proposed FE-CGFFT formulation was implemented last year (see Figs. 1 and 
2) but as can be expected, the rectangular mesh boundary does not always lead to the most 
efficient formulation, particularly when dealing with structures whose outer boundary is 
not rectangular. Because of this, during this year we developed and implemented a 
formulation which permits mesh termination at circular (see Fig. 3) boundaries for the 2D 
case with the corresponding boundary enclosure being a pillbox for the 3D case (see Fig. 

4). As before, these boundaries lead to convolutional integrals and do not therefore destroy 
the O(n) memory demand of the method. The FE-CGFFT formulation relating to circular 
(and ogival) boundary enclosures is described in the University of Michigan report 
025921- 11 -T (see also [13]) and results based on its implementation are shown in Figures 
5 and 6. Fig. 5 shows bistatic scattering patterns for a coated circular cylinder with a 
conductor radius of 3\, 0.05A. coating thickness and material properties £ \ = 3-j5 and (i r = 
1.5 - j0.5. The agreement with the series solution result is excellent. In Fig. 6 a 
backscatter pattern is shown for a V 2 x 1A. conducting ogive (see Fig. 3). In comparison 

with the moment method results, the agreement is again excellent. Additional results are 
given in Figure 7 for a missle-like shape scatterer. 


Pre- and Post-Processing Algorithms 

The availability of pre- and post-processing algorithms is crucial for the generation 
of the geometry and display of results in a graphical form. Generally, it is desirable that 
these tasks be done with a graphical user interface (GUI) and possibly in an X-window 
setting. Part of this year's effort was therefore devoted to the development of such 
algorithms and/or interfaces for the more sophisticated commercial pre- and post- 
processing packages. 
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For the most part, there exist commercial geometry, mesh generation and post- 
processing packages which are highly interactive and graphical. Nevertheless, there is 
always a need for a suitable interface or data interpreter between the commercial packages 
and the solution algorithm described in the previous section. The specific package 
interfaced with the computational algorithm was SDRC IDEAS@ and the selection of this 
was based on its availability on the U-M Network, its versatility, graphical user interface, 
and capability to generate meshes for 2D and 3D structures. Furthermore, a new version of 
SDRC IDEAS@, to be released soon, will support X- windows. IDEAS was developed 
for mechanical design purposes, but its geometry and finite element mesh generation 
modules are particularly suited for our needs. The geometry is defined graphically using 
the area (for 2D) or solids (for 3D) modeling capability provided by the module Geomod. 
Alternatively, the user may choose to enter the geometry in terms of individual points, 
surve segments (for 2D) or surfaces (for 3D). Once the geometry is entered, mesh areas or 
regions are specified and either a mapped mesh or free mesh can be generated. Further, 
individual nodes and/or elements may be inputted as desired using the CREATE command. 
Once the mesh is generated, two files are created, one containing the nodes and their 
corresponding coordinates, and another specifying the nodes of each element. These files 
are then read by an interpretor which creates a new input file compatible with the format 
required by the computational modules. 

Examples of two-dimensional meshes generated with SDRC IDEAS@ were shown 
in Figures 1, 2b and 7, and these are in a form suitable for the FE-CGFFT analysis. Some 
three-dimensional meshes are also displayed in Figure 8 for an ogive and missile-like 
structures. As seen, the 3D meshes are terminated at a cylindrical surface, tightly enclosing 
the scatterer which is the intended enclosure for the proposed FE-CGFFT method. A brief 
manual for geometry and mesh generation using SDRC IDEAS@ is currently being 
prepared. 

SDRC IDEAS@ is a rather sophisticated package and its use is certainly preferable 
for 3D modeling and mesh generation. For 2D mesh generation, though, it is possible to 
construct a non -proprietory package without much effort, and which is also simpler 
without a serious sacrifice in versatility. Clearly, the primary reason for resorting to such 
an algorithm is to permit mesh generation at sites not having a license for SDRC LDEAS@. 
The specific geometry and mesh generation package developed for this purpose is based on 
the algorithm described in [ 14]. Examples of free meshes generated by this package are 
displayed in Figures 2a and 3. The package is interactive/menu driven and can be readily 
used without much preparation. The mesh can be displayed in the Apollo screens or a 
postscript file may be generated for display on other workstations. At present. 
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visualization cannot be done in an X-window but this capability is planned for early next 
year. 

A variety of post-processing capabilities have also been exmployed for a graphical 
display of the output data. The output is either in the form of echowidth plots as a function 
of observation and/or incidence angle or in the form of gray level field maps. Color instead 
of gray level field maps can also be generated at those workstations which support this 
feature. Generally, all echowidth plots are generated and displayed using standard 
software, and each workstation provides its own selection. To generate and visualize the 
gray level and color field maps we employed the public domain package GRAFIC. An 
example of a gray level plot is shown in Figure 9. This is generated from a postscript file 
and can thus be displayed on other sponsor machines. 


3D Algorithm for Bodies of Revolution 

Before extending the presented 2D formulation to scattering by arbitrary 3D 
structures, it is instructive that we first consider its implementation for a restrictive class of 
3D bodies. In particular, during this year an algorithm was developed for scattering by 
inhomogeneous bodies of revolution. Because of the symmetry of this structure, it is only 
necessary to discretize it in a single plane slicing the structure as shown in Figure 10. A 
knowledge of the fields over this cross-section is then sufficient to generate the fields 
everywhere by employing a Fourier expansion in the azimuthal direction. Clearly, the 
discretization can be accomplished using a 2D mesh generation routine and this is the 
primary reasons for considering this class of structures. Also, the storage requirements are 
comparable to that of the 2D formulation although, as expected, the computational intensity 
is much greater. 

The mathematical details pertaining to the BOR formulation will be presented at 
sufficient detail in an upcoming technical report. Briefly, the method consists of the 
following steps 

1 ) A Fourier expansion is used to expand the fields in terms of those over a single 
cross-section of the BOR. 

2) The fields in the finite-element region are then formulated via the Coupled 
Azimuthal Potential (CAP) method as described in [15]. This results in a banded 
finite-element matrix in terms of the boundary fields. 
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3) The boundary fields are formulated via the usual Stratton-Chu equations which are 
then discretized via the boundary element method. As before, the boundary enclosure 
is chosen to yield convolutional integrals computed via the FFT. 

4) The finite-element and boundary-element systems are coupled via the boundary 
fields and solved via the CGFFT method maintaining an O(n) storage requirement, 
where n is the number of nodes over a single cross section of the BOR. 

Presently, a code has been written based on the proposed formulation and is in the 
final stages of the validation process. 


CONCLUSIONS 

: o far, we have formulated and implemented the FE-CGFFT method for a variety 
of 2D structures and we are now in the process of completing its implementation for BOR 
structures. The method was proposed because of its versatility, accuracy and low memory 
demand in comparison with other methodologies, and all of these attributes have been 
demonstrated in the testing and validation process. It is therefore a promising method for 
general 3D implementations to be considered in the following year. 


TRANSITIONS 

The validation of the 3D BOR formulation is expected to be completed by early Fall 
1990. We will then begin the development and implementation of the formulation for 
arbitrary 3D structures. This implementation is expected to be much more involved than 
those considered earlier and the same if true for the geometry and mesh generation. It is 
therefore, likely that the proposed 3D implementation may not be completed by the end of 
the 3rd year. Also, because of the need to generate suitable pre-processing and post- 
processing algorithms additional man-hours are required during the third year of this effort. 
Most likely, a practical user-oriented validated and benchmarked code will not be available 
until the fourth year. As part of this effort it would also be desirable to design and develop 
a graphical user interface (GUI) compatible with the X-window platforms. The GUI is 
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particularly necessary for the 3D analysis package. Otherwise, the user would be faced 
with a long list of subprograms whose interfacing would likely be cumbersome. 

The proposed FE-CGFFT formulation employs an exact boundary condition at the 
termination of the me' h. This eliminates a need to extend the mesh far from the scatterer 
leading to a substantial savings in storage requirements. However, this storage reduction 
and solution accuracy is achieved at the expense of computational complexity and intensity. 
In many cases, though, where accuracy is not of primary concern, one could resort to the 
use of non-exact (i.e. absorbing boundary conditions), for terminating the mesh. This 
leads to completely sparse matricies which can be solved n.ore efficiently using special 
purpose algorithms. In the future, it is therefore desirable to include this formulation as an 
option to the user. Also, a new class of boundary conditions are currently being 
investigated for terminating the mesh. 
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Figure 1 : E z and H z backscatter echowidth patterns for the illustrated 

perfectly conducting cylinder. 





COATED OGIVE 



Figure 2(b): Example of a rectangular mesh enclosing a coated ogive. 


FILE: ogive_leo_out 
CREATION DATE: 7 Aug 1990 
STRUCTURE: 0.5 x 1 X coated conductor 
ENCLOSURE: rectangle 

number of nodes: 1084 

number of elements: 1 936 

nodes on pec boundary: 140 

nodes on obs boundary: 92 

nodes on unknowns (E-pol): 944 
nodes on unknowns (H-pol): 1084 
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Figure 2(c): H-polarization backscatter pattern for the coated ogive in 

Fig. 2(b). 
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Figure 5: E z and H z bistatic echo width from a coated circular cylinder with 

a conductor radius of 3A, and coating thickness of 0.05X with 
material properites £ r = 5 - j5, p r = 1.5 - j0.5. 
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FILE: missile_out 
CREATION DATE: 1 Aug 1990 
STRUCTURE: missile shape 
ENCLOSURE: rectangle 

number of nodes: 636 

number of elements: 908 

nodes on pec boundary: 172 

nodes on obs boundary: 192 

nodes on unknowns (E-pol): 464 
nodes on unknowns (H-pol): 636 
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Missile Structure 



Figure 7: Backscatter patterns for a missile-like perfectly conducting and 

coated cylinder. 








Figure 8: 3D meshes for an ogive and missile-like structures. 








Figure 10: Body of revolution 


surrounded with a rectangular mesh. 
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